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Abstract. We computed ground-state energies of calcium isotopes from 42 Ca to 48 Ca by means of the 
Auxiliary Field Diffusion Monte Carlo (AFDMC) method. Calculations were performed by replacing the 
40 Ca core with a mean-field self consistent potential computed using Skyrme interaction. The energy of the 
external neutrons is calculated by projecting the ground-state from a wave function built with the single 
particle orbitals computed in the self consistent external potential. The shells considered were the IF7/2 
and the lF 5 / 2 - The Hamiltonian employed is semi-realistic and includes tensor, spin-orbit and three-body 
forces. While absolute binding energies are too deep if compared with experimental data, the differences 
between the energies for nearly all isotopes are in very good agreement with the experimental data. 

PACS. 21.10.-k Properties of nuclei; nuclear energy levels - 21.10.Dr Binding energies and masses - 
24.10.Cn Many-body theory - 21.30.Fe Forces in hadronic systems and effective interactions 
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1 Introduction 

The study of medium-heavy nuclei using a realistic Hamil- 
tonian is a very hard problem to attack because of the 
complexity of the nucleon-nucleon (NN) potentials and of 
the three-body forces (TNI) . The calculation of properties 
of a nuclear systems generally starts by considering an ef- 
fective Hamiltonian with a NN and TNI effective interac- 
tions describing the one-, two- and three-pion exchange 
between nucleons. Typically other effects are also included 
in the NN interaction, that is fitted to the experimental 
scattering data. 

Efficient methods have been developed to accurately 
solve the Schrodinger equation for few-nucleon bound states 
with a realistic NN interaction[ll[2] and with a three-body 
force[3,4J. However, nuclei which can be treated by these 
methods have very small mass numbers A(5J|5]. 

Projection Monte Carlo methods are the best candi- 
date to perform calculations with such potentials. In fact 
all the results obtained so far are in close agreement with 
experimental data. In turn, these procedures allowed for 
the development of more and more accurate model interac- 
tions that are now good enough to describe light nuclei:?] 
very accurately. 

All nucleons Green's Function Monte Carlo (GFMC) 
calculations are presently limited to systems with only 
12[8] nucleons or 149] neutrons. Medium-heavy nuclei are 
usually investigated with approximate many-body tech- 



Variational Monte Carlo 13], or other techniques [14lfT5j 
which are less demanding from the computational point 
of view, but which contain uncontrolled approximations. 

The study of neutron-rich isotopes is an interesting 
field of investigation because their structure is particu- 
larly relevant to constrain properties of the crust in neu- 
tron stars[16J. The impossibility to use accurate methods 
such the GFMC to study nuclei with intermediate mass 
values, in addition to the complexity of NN and TNI nu- 
clear forces imposes the need to introduce a model simpler 
than the full microscopic description of a nucleus includ- 
ing all degrees of freedom. The study of heavy isotopes 
of oxygen has already been performed within a scheme 
in which the closed-shell nucleus is substituted by an ex- 
ternal potential [17). In this paper we propose a similar 
analysis of heavy calcium isotopes in the IF7/2 and IF5/2 
shells based on the Auxiliary Field Diffusion Monte Carlo 
(AFDMC) method, that was already used in several works 
on pure neutron systems giving results of a quality com- 
parable with those of GFMC [IS]. This case presents some 
additional difficulties with respect to that of the oxygen 
isotopes. In fact the density of external neutrons is closer 
to that of 40 Ca, and therefore core polarization effects may 
become more important. The external neutrons, treated 
explicitly, interact with a realistic two and three body 
potential (in particular the Argonne AV8'[19 plus the 
Urbana IX ED]). Recently, the AFDMC was extended to 



niques like Hartree-Fock pUl I TT] . Correlated Basis Functions[l®bpton-neutron systems such nuclei 2 lj and nuclear matter [ 
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but using a Hamiltonian limited to include only the main 
contribution coming from the one-pion exchange. 

The aim of this paper is to verify whether the single- 
particle picture is still a good starting point for this kind of 
calculations and whether the AFDMC algorithm is able to 
account for the relevant quantum many-body correlations 
for neutron systems. In fact, for nuclei containing a num- 
ber of nucleons above the closed shell numbers of oxygen 
the effect of magic numbers might vanish, and the accu- 
racy of the shell model might be reduced. In addition, in 
this case, neutrons have a large angular momentum (I = 3) 
that might play an important role in the spin-orbit inter- 
action. It might be hard to correctly evaluate the contri- 
bution of the spin-orbit part to the total binding energy. 
In addition, if by replacing the core of 40 Ca with an ex- 
ternal well, we should find that the separation energies 
between different isotopes were in a good agreement with 
experimental data, then the contribution to energies sep- 
aration would come almost completely from NN and TNI 
used, and core polarization effects can be neglected. 

One could in principle include other shells having sim- 
ilar single particle energy, like the 2P 3 / 2 , in order to start 
from a more accurate trial wave function to project out 
the ground-state properties of the system. We will show 
that results obtained limiting the space to the IF shell 
are already very satisfactory. The inclusion of shells with 
principal quantum number different from 1 would intro- 
duce many additional technical difficulties [IT]. 



2 Hamiltonian and methods 



The ground-state energies of the calcium isotopes are com- 
puted starting from a non-relativistic Hamiltonian of the 
form: 

H = T + V 1 + V 2 + V 3 

= -E^ v * + E F - t ( r *)+E%' + E v *i", 

i i i<j i<j<k 

(i) 

where the particle index i and j labels only the N external 
neutrons. In our calculations we considered systems with 
N=2 to N=8. The one-body potential V ex t represents the 
40 Ca core; it has been obtain from Hartree-Fock calcula- 
tions using Skyrme forces. In fact the only input of the 
theory is the set of Skyrme parameters. Based on our pre- 
vious experience with oxygen isotopes we chose the set of 
parameters corresponding to the Skyrme I interaction [11]. 
However we expect that the choice of the self-consistent 
potential only influences the absolute total binding energy 
of each isotope. The energy difference between isotopes 
mostly depends on the choice of the NN and TNI interac- 
tions including correlations between the external neutrons. 
Our approximation corresponds to neglecting the effects 
of the interaction between external neutrons and the core. 



The two body interaction used belongs to the Urbana- 
Argonne family 23J; it is written as a sum of operators: 

l 

Vl(hj) = EE 

)OW(i,j). (2) 

In this work we truncated the sum to include only the 
first 8 operators with the parameters of the Argonne AV8' 
potential^!). This contains the usual operators: 

P=1 > 8 (l, j) = (1, CTi ■ (Tj, S ih Lii ■ Sij) X (1, Ti ■ Tj) , (3) 

where the operator Sij — icrt ■ rijCTj ■ fy — <r, • <jj is the 
tensor operator and = —ihrij x (Vj — Vj)/2 and 
Sij = h(cTi + <Tj)/2 are the relative angular momentum 
and the total spin for the pair ij. For neutrons T\-t j = 1, 
and we are left with an isoscalar potential. 

The AV8' potential used in this work is a simplified 
version of the more accurate Argonne AV18. It repro- 
duces very well the binding energy of nucleons for den- 
sities smaller then or equal to the equilibrium one. It has 
been used in calculations of light nuclei [2"U1|2"5] , symmet- 
ric nuclear matter [2"2~|26|, neutron matter [57], spin polar- 
ized neutron matter [2"B"] and neutron rich nuclei [T7]. The 
Hamiltonian contains also the Urbana IX (UIX) potential[20] . 

Calculations were performed using the AFDMC method 
|29j . This algorithm projects out the lowest-energy ground- 
state from a trial wave function ipT by a propagation in 
imaginary time r: 

iP(T)=e-( H - E ^ T , (4) 

and for sufficiently large r 

<p Q = lim t/)(r) , (5) 

r — >oo 

where 4>o is the lowest-energy component of tpT not orthog- 
onal to ipT- The evolution in imaginary-time is reached by 
solving the integral equation including importance sam- 
pling 

ifr(R)1>(R,r) = f dR'G(R,R',r)^§-^ T (R'WR\0), 
J Vt{R) 

(6) 

where G(R, R' , r) is the approximate Green's function of 
the system, that in the limit of small time-step is 

— 

/ m \ 2 m(R-R') 2 VIR)+V(R') . 

G{R,R! 1 At)=[ — ^ — e — ^^ e - 2 Ar . 
\2-kTi 2 At) 

(7) 

Then the integral of Eq. [6] must be iteratively solved until 
the convergence is reached. 

The AFDMC algorithm implements the usual diffusion 
process for the particle positions and samples the spin 
states of neutrons with a propagator written in terms of 
a Hubbard-Stratonovich transformation. This is done in 
order to reduce the quadratic spin dependence of nuclear 
Hamiltonians on the spin operators to an integral over 
auxiliary field variables, therefore averaging the sum over 
spin states. 
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As in standard DMC, AFDMC suffers of the fermion 
sign problem due to the antisymmetric character of the 
wave function; in our case the importance function is com- 
plex, and we constrain walkers to propagate within regions 
where the real part of the importance function has the 
same sign. A detailed description of the method and of 
the constrained path [30], used to control the fermion sign 
problem, can be found in[31j. 

2.1 Wave function 

The wave function used both as importance and projection 
function for the AFDMC algorithm has the following form: 

^(R,S)=Fj(R) D(R,S), (8) 

where R = (r±, . . . , rjy) and S = (s%, . . . , sjy). The spin 
assignments Si consist in giving the spinor components, 
namely 

Si = (jj^J = Ui\ T) + df| |) , 

where and di are complex numbers. 

The Jastrow correlation operator is given by 

Fj(R)=Hfj( rij ), (9) 

i<j 

while the antisymmetric part of the trial wave function is 

D(R,S) = A[^ aj {r i ,s i )] (10) 

that is the Slater determinant of one-body spin-space or- 
bitals: 

(f) a (r, a) = Rnj(r)Yi imi {9, 0)£ s , TOs (a) , (11) 

and a = {n, j, m,j}. 

The radial components R n _j(r) were obtained solv- 
ing the Hartree-Fock (HF) problem with the Skyrme I 
force [11]. The resulting self-consistent single particle po- 
tential is used in substitution of the closed 40 Ca core. The 
yielded radial functions are written in the j, rrij base. 

The angular components were chosen to be eigenfunc- 
tions of the total angular momentum operator J, in order 
to enforce the proper angular nodal structure and to re- 
duce the computational load. They were obtained by ap- 
plying projection operators of the required momentum on 
single-determinant wave functions, in the l—s representa- 
tion. Clebsch-Gordan coefficients were obtained automat- 
ically by projecting from single-particle wave function and 
constructing the transition matrix. 

In "Ca, n=42...48 isotopes, we assume that neutrons 
fill only the orbitals in the IF7/2 and I-F5/2 shell in order 
to build the ground-state of correct symmetry. The many- 
body states are obtained by coupling the single-particle 
angular momentum by constructing eigenstates of total 
angular momentum J = ji + ... +j'jv with N=2...8; for the 
states with an even number of neutrons, the ground-state 
has J = 0, while for odd neutron numbers, the ground- 
state has total angular momentum J = 7/2. These states 



are in general written in terms of a linear combination of 
Slater determinants, whose coefficients are determined by 
the symmetry of the state, and obtained by projection. 
Each determinant is evaluated at the current values of 
the positions and spin assignments of the nucleons in the 
walker \R,S). 

A sufficiently good representation of the ground-state 
of 46 Ca can be obtained by building a two- hole state which 
is complementary to 42 Ca. For the 43 Ca and 45 Ca com- 
pletely different trial wave function are needed; in fact 
the wave function for the first nucleus contains 9 deter- 
minants and the second 35. In principle we could use the 
43 Ca wave function as a three-hole state describing 45 Ca. 
However, AFDMC gives in this case an energy a few MeV 
too high. This is probably due to the approximation used 
to deal with the sign problem. The two wave functions are 
degenerate in the eigenspace of J 2 but they give different 
AFDMC energies, because of the different nodal surface. 

The Jastrow function fj has been taken as the scalar 
component of the Fermi Hypernetted Chain in the Single 
Operator Chain approximation (FHNC/SOC) correlation 
operator Fij which minimizes the energy per particle of 
neutron matter at density p=0.16 fm~ 3 [32]. The Jastrow 
part of the function in our case has the only role of reduc- 
ing the overlap of nucleons, therefore reducing the energy 
variance. Since it does not change the phase of the wave 
function, it does not influence the computed energy value 
in projections methods. For this reason the Jastrow func- 
tion has not been further optimized for our calculations. 

In order to have an accurate estimate of the ground- 
state energies of the isotopes we performed several sets 
of runs for different values of the imaginary time step 
and walker populations. The reported values are all ex- 
trapolated to At — > 0. Also the dependence of the result 
from the number of walkers used was investigated. Some 
calculations were repeated for 500 and 1000 walkers. For 
relatively long time-steps, the energy has rather large fluc- 
tuations when a smaller number of walkers is employed. 
On the other hand the average energy does not show a 
clear trend outside the statistical fluctuations. Therefore, 
we present all results obtained with 1000 walkers. 



3 Results 



isotope 


E ' AFDMC 


E eX p 


E( 4 "Ca) 


-25.25(7) 


-19.85 


E( 43 Ca) 


-33.9(4) 


-27.78 


E( 44 Ca) 


-46.3(1) 


-38.91 


E( 45 Ca) 


-52.6(4) 


-46.32 


E( 46 Ca) 


-62.9(3) 


-56.72 


E( 47 Ca) 


-70.2(7) 


-63.99 


E( 48 Ca) 


-80.3(8) 


-73.09 



Table 1. Ground-state energy computed with AFDMC. All 
the energies are given in MeV. Experimental values are referred 
to the ground-state energy of 40 Ca, taken from Ref. [33] . 
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Fig. 1. Outline of differences between energies of the isotope 
series studied. All the energies are expressed in MeV and all 
values are referred to the ground-state energy of 48 Ca from 
AFDMC and experiment respectively. 
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Fig. 2. Densities of external neutrons for all isotopes in the 
ground-state calculated with AFDMC and the Skyrme density 
of 40 Ca. 



In table [T] we report the AFDMC energies obtained for 
the isotopes series 42 Ca- 48 Ca, compared with the available 
experimental values [55]. 

As expected, the absolute binding energies are quite 
different from the experimental results, although the rel- 
ative discrepancy never exceeds 30%. This is a drawback 
of using of an external potential for including the effects 
of the filled core of the nucleus. The total binding energies 
are all overestimated. This reflects the absence of a correct 
description of the density of neutrons at the center of the 
drop, which is underestimated due to the absence of the 
S states, giving rise to an effective potential which is too 
deep for small distances of the neutrons from the center. 
Moreover we completely neglect core-polarization effects. 





Eafdmc 


E GX p 


E( 42 Ca)-E( 48 Ca) 


-55.1(8) 


-54.09 


E( 43 Ca)-E( 48 Ca) 


-46.4(8) 


-46.16 


E( 44 Ca)-E( 48 Ca) 


-34.0(8) 


-35.03 


E( 45 Ca)-E( 48 Ca) 


-27.7(8) 


-27.62 


E( 46 Ca)-E( 48 Ca) 


-17.4(8) 


-17.22 


E( 47 Ca)-E( 48 Ca) 


-10.1(8) 


-9.95 



Table 2. Ground-state energy differences (in MeV) computed 
by AFDMC among Ca isotopes considered in this work. 



Most of the information needed to understand the ef- 
fects of NN and TNI interaction in the external shell can 
be obtained looking at energy differences between the iso- 
topes considered. In fact, if the intra-shell interaction has 
a dominant effect, the gaps should not depend too much 
on the quality of the external well considered. 

In Tableland in Figure Q] we report the energy differ- 
ences for the isotope series considered, compared with the 
corresponding differences obtained from the experimen- 



tal results. As it can be seen, in this case the agreement 
between computed and experimental values is excellent. 
Small deviations are present only for the two cases of iso- 
topes 42 Ca and 44 Ca. 

In figureOwe report the AFDMC densities normalized 
to unity of the external neutrons for the isotopes consid- 
ered in this work. In the figure we also display the density 
of 40 Ca calculated with the Skyrme I force. As it can be 
seen the neutron's densities are all quite similar, and very 
small deviations are present. The external neutrons are 
very close to the core of 40 Ca and because of this one 
might expect that the interaction between the core with 
external neutrons cannot be a satisfactorily described by 
a one-body external potential. However, the results on the 
separation energies show that this effect only contributes 
at the single particle level. 

4 Conclusions 

We used the AFDMC method to study properties of off- 
shell neutrons in calcium isotopes including a realistic two- 
and three-body interaction. The approximation of consid- 
ering the core as a single particle interaction on external 
neutrons gives satisfactory results. In fact it does not affect 
the differences between energies of the isotopes considered, 
despite there is a systematic difference of about 6-7 MeV 
between our absolute energies with experimental values. 
This fact reveals that the physics of the external neutrons 
is dominated by the NN and TNI interactions. The impor- 
tance functions to project the ground-state with AFDMC 
for the isotopes considered have been obtained by restrict- 
ing the Hilbert space to the lF 7 / 2 and If 5/2 shells. The 
quality of our results suggests that other shells, which are 
technically harder to include do not need to be consid- 
ered and should introduce only very small deviations in 
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the results. The model we considered should be applied to 26. 
investigate the microscopical structure of heavier neutron- 
rich nuclei or of excited states. 27. 

We thank S. Fantoni and K. E. Schmidt for useful 
discussions. Calculations were partially performed on the 28. 
BEN cluster at ECT* in Trcnto, under a grant for su- 
percomputing projects, and partially on the HPC facility ^9. 
"WIGLAF" of the Department of Physics, University of 
Trento. 30 - 
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